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Abstract 

A detailed study is presented of the search for Higgs bosons in the b-decay channel 
in the central exclusive production process at the LHC. We present results for proton 
tagging detectors at both 220m and 420m around ATLAS or CMS. We consider two 
benchmark scenarios; a Standard Model (SM) Higgs boson and the m™'^^ scenario of 
the Minimal Supersymmetric Standard Model (MSSM). Detector acceptance, smearing 
and event trigger strategies are considered. We find that the SM Higgs will be chal- 
lenging to observe in the b-jet channel without improvements to the currently proposed 
experimental configuration, but a neutral scalar MSSM Higgs Boson could be observable 
in the b-jet channel with a significance of 3a or greater within three years of data taking 
at all luminosities between 2 x 10^^ cm~^ s~^ and 10^^ cm^^ s~^, and at 5cj or greater 
after three years in certain scenarios. 
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1 Introduction 



The physics potential of forward proton tagging at the LHC has attracted a great deal of 
attention in recent years [Il[2l[3llll[5l[6l[71[9l[ini[IIl[l2l[l3l[Il]. A main focus of interest 
is the central exclusive production (CEP) process, pp ^ p + (p + p, in which the protons 
remain intact and the central system 4> is separated from the outgoing protons by a large 
rapidity gap. To a very good approximation, (p is constrained to be in a colour singlet, 
Jz = 0, state. Observation of any particle, such as a Standard Model Higgs boson, in 
the central exclusive channel would therefore provide a direct observation of its quantum 
numbers. Furthermore, by detecting the outgoing protons and measuring their energy 
loss accurately, it is possible to measure the mass of the centrally produced particle 
regardless of its decay products [15j. Because of these unique properties, it has been 
proposed that forward proton detectors should be installed either side of the interaction 
points of ATLAS and/or CMS. The FP420 collaboration has proposed to install detectors 
in the region 420m from the interaction points \16\ [T7] . These detectors would allow the 
detection of central systems in the approximate mass range 70 GeV < M^ < 150 GeV. 
Proposals also exist to upgrade the capabilities of ATLAS and CMS to detect protons 
in the 220m region [18^ 119] . These detectors, when used in conjunction with 420m 
detectors, would extend the accessible mass range well beyond 150 GeV. 

In this paper, we focus on the central exclusive production of the Standard Model 
(SM) Higgs boson and a Supersymmetric (MSSM) Higgs boson, with m^ ~ 120 GeV. 
The CEP process is shown schematically in figure [TJ For this mass region, the dominant 
decay channel of the Higgs boson is to bb, which is very difficult to observe in conventional 
Higgs searches at the LHC because of the large QCD background. This is not the case 
in central exclusive production due to the Jz = selection rule, which suppresses the 
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leading order central exclusive bb background by a factor of where M is the mass 
of the bb di-jet system. As we shall see, this renders the bb decay channel observable at 
the LHC in certain scenarios if appropriate proton tagging detectors are installed. 

The structure of the paper is as follows. Firstly, we give a brief overview of the 
proposed forward detector upgrades at 220m and 420m at ATLAS and CMS, including 
a simulation of the acceptance of the detectors. We then discuss the predicted signal 
cross sections and survey the background processes. Taking a 120 GeV Standard Model 
Higgs boson as the benchmark, we perform a simulated analysis including an estimation 
of the detector acceptance and smearing effects and Level 1 trigger strategies. The 
analysis is then extended to the MSSM for a particular choice of parameters. 

2 Forward detectors at the LHC 

There are two properties of the proposed forward detector systems that are critical to this 
analysis; the acceptance of the detectors in the mass range of interest and the ability of 
the forward detectors to correctly associate the detected outgoing protons with a Higgs 
boson candidate event measured in the central detector. This matching is critical at 
high luminosities, where the large number of proton-proton collisions per bunch crossing 
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Figure 1: Central exclusive production of a Higgs boson. 



(often referred to as pile-up) leads to a high probability that forward protons from single 
diffractive or double pomeron (DPE) collisions not associated with a Higgs candidate 
event will enter the forward detectors during the same bunch crossing. The proposed 
forward detectors aim to associate particular protons with the Higgs candidate event by 
making a measurement of the outgoing proton time-of-flight (TOF) from the interaction 
vertex to the detectors. The difference in the arrival times of the protons on opposite 
sides of the central detector, ATOF, allows a vertex measurement to be made under the 
assumption that the detected pair of protons originate from a single hard interaction. 
This vertex can then be matched with the vertex of a candidate Higgs event reconstructed 
using the central detector alone. The current design goal of the forward detectors is to 
achieve a lOps accuracy in the TOF measurement [16] . which translates into a vertex 
measurement accurate to 2.1mm. The use of fast-timing measurements is discussed in 
more detail in section 14.31 It has been suggested that the central detector could also 
be used provide a third timing measurement [20], which could allow for an improved 
rejection of pile-up events. We discuss the effect of this possibility in section [6l 

The acceptance of the forward detectors is governed by the distance of the active 
edge of the detector from the beam, which determines the smallest measurable energy 
loss of the outgoing protons, and the aperture of the LHC beam elements between the 
interaction point and the forward detectors; protons that lose too much momentum will 
be absorbed by beam elements, imposing an upper limit on the measurable momentum 
loss of the protons. 

The distance of the active edge of the detector from the beam depends primarily on 
the beam size at each detector location. Previous estimates [H] have assumed that the 
closest distance, D, is given by 

D w lOcJf, + 0.5mm (1) 

where cjfe is the gaussian beam size at the detector location and 0.5 mm is a constant 
term that accounts for the distance from the sensitive edge of the detector to the bottom 
edge of the window. The beam size at is approximately 250 ^m at 420m and 100 fim 
at 220m, leading to a distance of closest approach of 3 mm for detectors at 420m and 
1.5 mm for detectors at 220m. It is likely that the detectors will begin operation at 
a larger distance from the beam, at least until the detectors and machine background 
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Figure 2: The acceptance of the forward detectors as a function of central system mass 
for both protons tagged at 420m (a) or an asymmetric tag of one proton at 220m and 
one at 420m (b). The dependence of the acceptance on the distance of each detector 
from the beam is shown by the different curves. 



conditions are well understood |17l I22j . Figure [2] (a) shows the acceptance for events in 
which both outgoing protons are detected at 420m around IPl (ATLAS), as a function 
of the mass of the central system for different detector distances from the beam. The 
protons are generated using the ExHuME Monte Carlo [23] and tracked through the 
LHC beam lattice using the FPTrack program with version 6.500 of the LHC optics [24J. 
The central system mass, M, is calculated from the forward proton momenta, 

M2 ~ 6 ^25, (2) 

where the are the fractional momentum losses of the protons and ^/s is the centre-of- 
mass energy of the collision. The acceptance for a 120 GeV Higgs boson is independent 
of the distance of approach of the detectors from the beam up to approximately 7mm in 
the case where both protons are detected at 420m. This is consistent with the findings 

of [m. 

The situation is very different for events in which one proton is tagged at 220m and 
the other at 420m, as shown in figure[2](b). In this case, the acceptance is increased when 
either detector is moved closer to the beam. For a 120 GeV Higgs boson, with 420m 
detectors at 5mm and 220m detectors at 2mm, the acceptance is 28% if both protons 
are tagged at 420m (symmetric tag) with an additional 16% acceptance if one proton is 
tagged at 220m and one at 420m (asymmetric). Moving the 420m detectors inwards to 
3mm and the 220m detectors to 1.5mm increases the asymmetric acceptance by up to 
a factor of three, as shown in figure [2] (b). Figures [2] (a) and (b) also demonstrate the 
increasing importance of 220m detectors as the mass of the central system increases. At 
IP5 (CMS) the symmetric acceptance is identical to that at IPl. For the asymmetric 
tags, however, the acceptance is worse by a factor of 3 across the mass range of 
interest. This is caused by the horizontal (rather than vertical) plane of the crossing 
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angle of the beams at IP5 |17j . In this paper we concentrate on the measurement around 
IPl (ATLAS). 

3 Central exclusive H ^ bb production 

Central exclusive signal and background events are simulated with parton showering 
and hadronisation effects using ExHuME vl.3.4 [23]. ExHuME contains a direct imple- 
mentation of the Khoze, Martin and Ryskin (KMR) calculation of the central exclusive 
production process [3l[25]. The cross section for the CEP of a Standard Model Higgs 
boson decaying to bb as a function of the Higgs mass is shown in figure O The rapid 
decrease in cross section at m/j ~ 160 GeV occurs because the primary decay channel 
changes from bb to For masses above 140 GeV, it is expected that the Higgs 

boson should be observed in the M^Ty*-*-* channel with a luminosity of 30 fb^^ using 
forward proton tagging [26j . 

The primary uncertainties in the predicted cross section come from two sources; the 
parton distribution functions (pdf) and the soft survival (often termed rapidity gap sur- 
vival) factor. The CEP cross sections are relatively sensitive to the pdfs because the 
derivative of the gluon density enters to the fourth power. Figure [3] shows the cross sec- 
tion prediction for three different pdf choices, CTEQ6M, MRST2002NLO and CTEQ6L. 
The cross section for the CEP of a Standard Model Higgs boson of mass mh = 120 GeV 
decaying into b-quarks varies from 1.86 fb (MRST2002NLO) to 7.38 fb(CTEQ6Ll). The 
spread of a factor of ~ 4 is consistent with the findings of [9]. For the purposes of this 
analysis, we chose CTEQ6M as our default pdf as it lies between the two extremes. There 
is some justification for choosing an NLO pdf because the KMR calculation contains a 
NLO QCD K-factor (1.5) for SM Higgs Boson production. 

The soft survival factor, S^, is the probability that there are no additional hard 
scatters in a single pp collision and is expected to vary from process to process. For 
CEP processes, we take the ExHuME default of 5^ = 0.03 at LHC energies. Until very 
recently, there was a consensus that should be between 2.5% and 4.0% for the CEP 
of a Higgs boson at the LHC |27j . Two very recent studies have predicted a lower value 
|28 t 129 ^ 130] and it remains to be seen whether a new theoretical consensus can be reached. 
In any case, S'^ will be measurable in early LHC data. 

The MSSM contains three neutral Higgs bosons - two scalar and one pseudo-scalar. 
The Jz = 0, parity-even selection rules strongly suppress CEP of the pseudo-scalar. 
This can be advantageous in areas of MSSM parameter space where the pseudo-scalar is 
almost degenerate in mass to one (or both) of the scalar Higgs bosons, since CEP would 
provide a clean and complimentary measurement of the mass of the scalar only [3Tj , and 
allow nearly-degenerate Higgs Bosons to be distinguished [11]. Furthermore, at large 
values of tan/J, the cross section for the CEP of the scalar Higgs bosons can be strongly 
enhanced relative to the CEP of a SM Higgs boson of the same mass ^32j. We choose 
a point in parameter space defined by the m™'^^ scenario [33] with ruA = 120 GeV and 
tan/3 = 40, resulting in the scalar Higgs boson having a mass of 119.5 GeV and a decay 
width of 3.3 GeV. With this choice of parameters, the lightest scalar Higgs boson has an 
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Figure 3: The cross section for the central exclusive production of a Standard Model 
Higgs Boson decaying to b-quarks as a function of nih, for three different proton parton 
distribution functions. 

enhanced CEP cross section and is almost degenerate in mass with the pseudo-scalar. 
The CEP cross section for the lightest Higgs Boson decaying to b-quarks in this scenario 
is predicted to be 20.5 fb. The uncertainty on this prediction is the same as that for the 
SM Higgs boson. 

4 Backgrounds to the H ^ bb channel 

The backgrounds to central exclusive Higgs production can be broken down into three 
categories; central exclusive di-jet production, double pomeron exchange and overlap. 

4.1 Central exclusive di-jet backgrounds 

The most difficult backgrounds to deal with are the central exclusive (CEP) di-jet back- 
grounds. Central exclusive bb production is suppressed at leading order by the Jz = 
selection rule, but is still present and forms an irreducible continuum background beneath 
the Higgs mass peak. Central exclusive glue-glue production has a much larger cross 
section than bb production since it is not suppressed. It contributes to the background 
when the gluon jets are mis-identified as b-jets (1.3% probability for each mis-tag at AT- 
LAS - see section I5.3.2p . Central exclusive cc events do not contribute to the background 
for two reasons. Firstly, cc production is suppressed with respect to the bb process by 
^ . Secondly, the cc background is further suppressed by a factor of 36 relative to the bb 
cross section assuming that the probability of mis-identifying a c-quark is approximately 
0.1 for a b-tag efficiency of 0.6 [3^ . 

Higher order CEP di-jet backgrounds such as the 3-jet final state bbg process have 
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Figure 4: Production of a bb system via double pomeron exchange (DPE). 

been studied in [35j . It is expected that these backgrounds will be lower than the LO gg 
and bb backgrounds after experimental cuts, but they should be implemented into the 
ExHuME Monte Carlo and a full study performed before definitive conclusions can be 
drawn. 

4.2 Double Pomeron backgrounds 

Double pomeron exchange (DPE) is defined as the process pp ^ p + X + p, where X is 
a central system produced by pomeron-pomeron fusion (PP). The pomeron is assigned 
a partonic structure and so there are always pomeron remnants accompanying the hard 
scatter as shown in figure [H The relevant background processes are separated into bb 
and jj, where j represents light-quark and gluons jets. 

The DPE events are simulated using the POMWIG v2.0 event generator [HB], which 
implements the diffractive parton distribution functions measured by the HI Collabo- 
ration \37\ I38j . We use HI 2006 fit B, although we discuss the effect of using different 
diffractive pdfs. We choose the POMWIG option to treat the valence partons in the 
pomeron as gluons. POMWIG is normalised to the HI data, and therefore does not 
account for the soft survival factor. For hard double pomeron events, we take 52 = 0.03. 
POMWIG is also capable of generating scatters involving sub-leading (non-diffractive) 
exchanges (IRIR). We do not include this contribution because the cross section for 
IRIR ^ bb events is ~ 0.034 fb for ^1,2 < 0.05, 80GeV < M < 160 GeV arid 6-quark 
Pt ^ 35 GeV. This is negligible compared to the cross section for PP bb, which is 
~65 fb in the same kinematic region. 

4.3 Overlap backgrounds 

Overlap events are defined as a coincidence between an event that produces a central 
system of interest, X, and one or more diffractive events (single diffractive pp — > pX 
or DPE) which produce protons in the acceptance range of a forward detector in the 
same bunch crossing. On average there will be 3.5 interactions per bunch crossing at 
low instantaneous luminosity (10^^ cm~^ s~^) and 35 interactions at high instantaneous 
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luminosity (10^^ cm~^ s~^). We investigate three types of overlap event; [p][X][p], [pp][X] 
and [pX] [p] . The square brackets specify the interaction to which a part of the overlap 
event belongs - in this notation both the CEP and DPE events would be [pXp]. 
The cross section, aoiap, for the overlap background may be estimated by 



olap — cr [X] 



N=3 ' N=2 



\N - 
N=2 

where a^x] is the inclusive {pp — > jets) di-jet cross section, A is the average number of pp 
interactions per bunch crossing and N is the actual number of interactions in a specific 
bunch crossing. Because the actual number of interactions is not fixed, we sum over 
all possible numbers and weight each configuration by a Poisson distribution. In the 
first term, P2[p] (n) is the probability that, given n interactions, there are at least two 
events that produce a forward proton (one on each side of the interaction point) by any 
mechanism. This is dominated by soft single diffractive events pp pX. P2[p] (™) is 
given by a trinomial distribution, i.e. 

n r+q—1 , 
^ — ^ X — ^ n' / \^'/ \ Q / \n—r—q 

pmi-)= E E (/a) (/bO (i-/;i-/bi) w 

r+q=2 q=l ^ ' 'i'r ■'i- 

where, for example, /j^j is the fraction of events at the LHC that produce a forward 
proton in the +z direction and within the forward detector acceptance. 

In the second term, P[pp](n) is defined as the probability that there is at least one 
event that contains an outgoing proton on each side of the interaction point within the 
acceptance of the forward detectors (dominated by the soft double pomeron process 
pp PpX). P[pp](n) is a binomial distribution that utilises the event fraction /[ppj. Note 
that double pomeron events can also contribute to the first term in equation (Sj in the 
case where only one of the outgoing protons falls within the detector acceptance. 

The third term deals with a two- fold coincidence between a single diffractive di-jet 
event [pp — > pX), which produces a proton within the forward detector acceptance AND 
a hard central diffractive system X that mimics the signal, and an overlap event that 
produces a proton on the opposite side. <T[pX] is the total single diffractive di-jet cross 
section [pp — > pX)^ i.e the outgoing proton can be on either side of the interaction point. 

(n) is defined as the probability that there is at least one event with a forward proton 
in the forward detector acceptance on the opposite side of the IP to the single diffractive 
proton from the hard event. This is defined in a similar way to P[pp](n), but using the 
event fraction /[p]. 

For small proton momentum losses, the cross section for events at the LHC that 
contain a forward proton is dominated by the single diffractive cross section, usd [39j; 
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(a) (b) 

Figure 5: Figure (a) shows the di-jet [p][X][p] background cross section for both protons 
tagged at 420m. The di-jet cross section is calculated using HERWIG + JIMMY for 
partons with Et > 40GeV. Figure (b) shows the cross sections for each of the three 
types of overlap background. 



where ctt is the total cross section at the LHC, Sgjj is the soft-survival factor for soft 
single diffraction (0.087), ap(t) is the pomeron trajectory, t is the squared 4-momentum 
transfer at the proton vertex, gjyit) is the pomeron nucleon coupling and g^^p is the triple 
pomeron vertex. The value of /[p] is calculated by integrating equation [5] using standard 
Monte Carlo techniques, for a specific acceptance range in ^ and t. The cross section for 
the [p][X][p] di-jet events, for parton pr > 40GeV and 0.002 < ^ < 0.02, -5 < t < is 
shown in figure [5] (a). The cross section increases by two orders of magnitude from low 
to high luminosity. 

At larger values of ^, there will be an additional contribution from non-diffractive (i.e. 
reggeon exchange) events. This contribution is estimated using the PYTHIA [lO] and 
PHOJET [¥T] event generators. The predictions for the diffr active and non-diffractive 
contributions to /[p] are compared to equation [5] in table [H As expected, the diffractive 
contribution dominates at small The non-diffractive contribution becomes increasingly 
important at higher ^, but remains smaller than the diffractive contribution for ^ < 0.05. 
PHOJET predicts a higher fraction of non-diffractive events than PYTHIA and also 
predicts that a large fraction of forward protons are due to DPE events. In this analysis 
we use the prediction for single diffraction given by equation [5] for ^ < 0.02 (which would 
give hits only at 420m), and take the non-diffractive contribution to be negligible. For 
0.02 < ^ < 0.2, we add the PYTHIA prediction for the non-diffractive component, since 
this agrees more closely with theoretical expectations that the non-diffractive fraction of 
pp collisions at the LHC in this kinematic range should be between 1.0% and 1.7% [42j . 

The fraction of soft double pomeron events (DPE) within the detector acceptance, 
/[pp], is estimated using the PHOJET event generator. PHOJET predicts that the 
fraction of DPE events that have both forward protons in the range 0.002 < ^ < 0.02 
is 6.3 X 10~^. The fraction of DPE events that have one proton with 0.002 < ^ < 0.02 
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^ ran 


ge 






0.002 - 0.02 


0.02 - 0.05 


0.02 - 0.1 


0.02 - 0.2 


Equation 5 (fso) 


0.0112 


0.0040 


0.0070 


0.0098 


PYTHIA (fso) 


0.0104 


0.0045 


0.0081 


0.0112 


ifND) 


0.0002 


0.0016 


0.0043 


0.0124 


PHOJET ifsD) 


0.0069 


0.0031 


0.0055 


0.0081 


ifsD + IdPe) 


0.0097 


0.0045 


0.0081 


0.0118 


Und) 


0.0018 


0.0025 


0.0059 


0.0192 



Table 1: The fraction of events at the LHC that produce a forward proton on one side 
of the interaction point in a specific kinematic range. The PYTHIA and PHOJET event 
generators are compared to the single diffractive cross section given in equation [5l SD 
labels the outgoing proton from single diffractive scatters and ND labels the protons 
produced from non-diffractive scatters. DPE labels double pomeron exchange events. 



and one with 0.02 < ^ < 0.2 is 1.4 x 10^^. This allows us to construct the [pp][X] 
overlap backgrounds. The cross sections for the three types of di-jet overlap backgrounds, 
[p][X][p], [pX][p] and [pp][X] are shown in figure [5] (b) as a function of instantaneous 
luminosity. 

Monte Carlo overlap events used for the simulated analysis presented in section [5] are 
constructed as follows: The HERWIG Monte Carlo [13] is used to generate di-jet events, 
interfaced to JIMMY [3^ to simulate secondary scatters between spectator partons in 
the interacting protons. Single diffractive di-jet events are generated using POMWIG. 
For small diffractive protons are generated according to the distribution in ^ and 
t given by equation and added to the event record. For ^ > 0.02, the protons are 
selected as diffractive or non-diffractive on an event-by-event basis, according to the 
cross sections presented in section 14.31 We model the non-diffractive proton kinematics 
using a Regge flux factor, given by 

(fa e^^* 
where a^it) = 0.5 + 0.3t and bu = 1.6 GeV'^ [371 [38] • 

5 Simulated measurement of central exclusive Higgs Boson 
production in the b-jet channel 

Having generated the signal and background events, we now investigate a possible exper- 
imental strategy for detecting central exclusive Higgs production in the bb decay channel. 
The analysis presented here is specific to the ATLAS detector, in that we use the ATLAS 
parameters for detector smearing effects and b-tagging efficiencies. We do not expect 
the conclusions to be materially different for CMS for 420m detectors alone. However, as 
discussed in section [2] the asymmetric events, with one proton detected at 220m and one 
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at 420m, have a lower acceptance at CMS due to the different LHC optics. Ah figures 
are presented for the symmetric tag case, i.e. a proton hit within the acceptance of a 
420m detector on each side of the central detector 0. We present numerical results for 
both symmetric and asymmetric (420m + 220m) tags in section [6l 

5.1 Triggering on CEP b-jet events 

Proton tagging detectors at 420m are too far away to be included in the level one (LI) 
trigger decision at ATLAS or CMS in normal running and with the standard trigger 
hardware. It will therefore be necessary to design a trigger strategy that can keep the 
symmetric 420m + 420m events. Such a strategy must rely only on information from 
the central detectors. Information from 220m detectors can be included at level one, 
leading to a higher trigger efficiency for asymmetric 220m + 420m events. 

At ATLAS, the LI trigger threshold for jets is foreseen to be in the region of 180 
GeV and 290 GeV at low and high luminosity respectively. The reason for the high Et 
thresholds is that the QCD jet event rate is large and there is little additional rejection 
available at level 2 (L2), which can output a maximum rate of ~ 2 kHz at ATLAS. 
It is possible, however, to incorporate information from 420m detectors into L2. The 
requirement that there are opposite-side forward proton hits at 420m reduces the event 
rate significantly, as shown in figure [5] (a). Furthermore, if timing detectors are used to 
fix the vertex position, there is an additional rejection of a factor of ~ 20 as shown in 
section 15.5.11 The event rate at the output of L2 can therefore be reduced by requiring 
two in-time 420m proton hits at L2 by factors of 21500, 570 and 140 for luminosities 
of lO'^^, 5x10^^^ and 10^^ cm~^ s~^ respectively. A strategy could be to allow CEP a 
fixed (pre-scaled if necessary) input rate at LI for di-jet events with Ej- > 40 GeV. This 
high LI rate would then be reduced to an acceptable L2 output rate by proton hit and 
timing information in the 420m detectors. Further reductions in rate would be possible 
by cutting on the final state jet variables defined in section [221 

We present results as a function of maximum allowed input LI jet trigger rate. The 
maximum LI input rate we consider is 25kHz, which we label J25. This trigger would 
be un-prescaled for 40 GeV di-jets at 10^'^ cm~^ s~^, and be pre-scaled by a factor of 10 
at 10^^ cm~-^ s~^. Given the rejection factors above from requiring two in-time proton 
hits at L2, the maximum L2 output rate at 10^^ cm~^ s~^ would be 150Hz. This can 
be further reduced to a few Hertz at L2 output by applying cuts on the final state jet 
topology. We also present results for a lower LI input-rate of lOkHz (JIO). 

In addition to the jet trigger, low transverse momentum muon triggers will also save 
b-jet events. At low luminosity, a single un-prescaled muon trigger with pt ^ 6 GeV 
is possible. At high luminosity the pT threshold would be increased to reduce the rate. 
The rates for the muon triggers could be reduced significantly by requiring a coincidence 
with a jet with Et > 40 GeV. We use the notation MU6 for a 6 GeV muon trigger 
and MUlO for a 10 GeV muon trigger. We assume that no pre-scale is applied, which 

^We apply a cut of 80 < AM < 160 GeV on the mass M of the central system as measured in the 
forward proton detectors to restrict the size of the Monte Carlo samples before making the figures. 
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may mean that the jet requirement is included in the trigger at high luminosity. MUG 
will save 11% of b-jet events. If the threshold were increased to 12 GeV (MU12), 4% of 
events will be saved. 

At low luminosity, a rapidity gap requirement between the outgoing protons and the 
central system in addition to two central jets with Et > 40 GeV will have an acceptable 
trigger rate. This trigger is self-prescaling as luminosity increases, because pile-up events 
will destroy the rapidity gaps. Using a Poisson distribution, we evaluate the probability 
that there is only one inelastic scatter in the bunch crossing, i.e the event of interest. 
This rapidity gap trigger provides a luminosity-dependent trigger efficiency of ~ 17% at 
10^3 cm-2 s"^ and ~ 2% at 2 x lO^^ ^m-^ s'^. 

As mentioned above, the 220m detectors can be included in an LI trigger decision. 
For instantaneous luminosities up to 2-3 xlO'^^ cm^^ s^^, it is expected that a trigger 
consisting of one proton hit at 220 m and two 40 GeV central jets requires no pre-scale 
to have a LI trigger rate of IkHz [45]. If, as for the di-jet triggers, a relatively large LI 
input trigger rate is allowed which can be reduced at L2 by 420m hit information, then it 
should be possible to achieve a trigger efficiency close to 100% for the asymmetric events 
at high luminosity, although further studies will be required to confirm this expectation. 

The efficiencies of the above trigger strategies and the effect on the significance of 
the measurement are presented in section [62J 



5.2 Selection of exclusive di-jet events in the central detector 

Experimentally, the definition of a central exclusive di-jet event, i.e. an event with no 
activity outside the jets, is somewhat arbitrary. The CDF collaboration have searched 
for a central exclusive di-jet signal [1^ using the di-jet mass fraction, Rjj, where 

Rn - -Jf^ (7) 

where Mjj is the mass of the di-jets and M is the mass of the central system given by 
equation [2]EI- Exclusive events might be expected to have a di-jet mass fraction close 
to 1, since there should be no activity outside the jets. However, parton showering, 
hadronisation and out-of-cone effects (CDF use a cone-based jet algorithm) smear the 
measurement, so that even for true exclusive di-jet events, Rjj can be significantly less 
than 1 [i7| . 

An alternative definition of the di-jet mass fraction was proposed in [IS] which aims 
to reduce the effects of hard final state radiation. In a central exclusive event, the highest 
transverse energy jet is the one that has been least affected by final state radiation. The 
di-jet mass fraction Rj is defined as 

Rj = ^cosh{r]i-y), (8) 



^CDF do not have proton taggers on each side of the experiment, and therefore estimate and ^2 



using the central detector. Details can be found in |4Bj 
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where and rji are the transverse energy and pseudo-rapidity of the leading jet in the 
event. The pseudo-rapidity of the central system, y, is obtained from the measurement 
of the forward protons; 



The pseudo-rapidity of the central system can also be obtained from the average pseudo- 
rapidity of the two jets. This means that the variable 



should be approximately zero for a central exclusive event. 
5.3 Jet algorithms 

In this section we investigate the dependence of the jet-based variables Rjj, Rj and Ay 
on the choice of jet algorithm. We use the mid-point cone algorithm [49j with varying 
cone radius and the Kt algorithm [50] with different values of the R-parameter. Events 
are accepted if the highest transverse energy jet has Ej- > 40 GeV. 

5.3.1 Generator level 

The Rjj distributions are shown in figuresE](a) andEKb) for the cone and Kj- algorithms 
respectively. As expected, the distributions peak at Rjj ~ 1.0 with a tail that extends 
to lower values due to parton showering and hadronisation effects. As the cone radius, 
R, is decreased, the distribution is smeared to lower values of Rjj because for smaller 
cone radii, particles are more likely to be reconstructed outside of the jets and the 
reconstructed mass of the di-jet system is correspondingly lower. A similar result is 
found for the Kt algorithm when reducing the R-parameter, RKrp. 

The dependence on the cone radius and R-parameter is less marked for the Rj vari- 
able. The distributions are still peaked at Rj ~ 1.0, but changing the cone radius or 
R-parameter has only a small effect on the distribution as shown in figures (c) and E] 
(d). This is expected because the leading jet is less affected by final state radiation, and 
should therefore be less affected by changes in jet definition. 

The Ay distributions are shown in figures [6] (e) and [6] (f) for the cone and Kt algo- 
rithms respectively. The distributions are strongly peaked at zero for both jet algorithms 
and all the choices of cone radius and R-parameter result in the majority of events being 
reconstructed in the region Ay < 0.06. We conclude that the Ay and Rj variables are 
relatively insensitive to the choice of jet algorithm before detector effects are taken into 
account. 

5.3.2 Detector effects 

Figure [7] (a) and (b) show the effects of detector smearing on the Rj and Ay distributions 
for CEP Higgs events (i.e. signal only) for a cone radius of 0.4. The particle energy. 




(9) 




(10) 
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Ay Ay 



(e) (f) 

Figure 6: The dependence of the Rjj, Rj and Ay variables on the jet algorithm used to 
reconstruct the jets. Figures (a), (c) and (e) show the dependence of these variables on 
the radius used in the cone algorithm. Figures (b), (d) and (f) show the dependence on 
the Kt di-jet. 
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Figure 7: The effect of detector resolution on the Rj and Ay variables of the recon- 
structed signal events. The resolution of the ATLAS calorimeters dominates the smear- 
ing of the Rj variable (figure (a)) and the resolution of FP420 only has a small effect. 
Figure (b) shows that the Ay variable is only slightly affected by the detector resolution. 



pseudo-rapidity and azimuth are smeared by the resolution of the ATLAS electromag- 
netic and hadronic calorimeters [51] • The momentum of the leading protons measured 
by the forward detectors is smeared by the resolution of the 420m detectors [12] . 

The Rj distribution is smeared around the peak value. The effect of the ATLAS 
central detector resolution on the Et measurement of the leading jet is the dominant 
effect, with the smearing of the mass measurement from the 420m detectors contributing 
only a small amount. For the Ay distribution, the detector effects are minimal. 

In addition to smearing, the 6-tagging efficiency of the ATLAS detector must be 
taken into account. The single 6-tagging efficiency at ATLAS is 60% for a 6-jet, with a 
mis-tag rate of 1.3% for jets originating from gluons [3l]. Therefore the efficiency for a 
system consisting of two 6-jets being experimentally identified as a two 6-jet system is 
0.36. The mis-tag rate for a gluon and b di-jet is 0.0078 and the mis-tag rate for a two 
gluon di-jet is 1.69xl0~^. The c-quark mis-tag rate is 0.1. 

5.4 Identification of diflFractive events 

Diffractive events have historically been identified by searching for a rapidity gap between 
the central system and the outgoing proton. This approach has a high efficiency at the 
LHC only at very low luminosities, (around lO'^'^ cm~^ At higher luminosities, 

particles from pile-up events will in most cases destroy the rapidity gap. As stated in 
section [571] the probability for no inelastic pile-up events is ~17% at an instantaneous 
luminosity of 10^^ cm"^ s~^, reducing to ~2% at 2 x 10^^ cm~^ s~^. 

Because of this, we use charged track multiplicity cuts to identify CEP candidate 
events. The primary hard scatter in a non-diffractive overlap event is characterised by 
a central system that is colour connected to the proton remnants. There will therefore 
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be on average more charged particles associated with the di-jet vertex in an inelastic 
collision than in a CEP collision. Furthermore, the underlying event, which cannot be 
present in a CEP event, will be an additional source of charged particles at the vertex. 
At ATLAS, the inner detector has excellent charged track vertex reconstruction, which 
allows charged tracks to be identified with a specific di-jet vertex. We use two variables 
to identify CEP candidate events using the charged particles that are associated with 
the di-jet vertex: 

• The number of charged particles, Nc-, that are not contained within the two highest 
Et jets. Nc is a measure of the particle multiplicity in the event that is not 
associated with the hard scatter. It is of course dependent on the jet algorithm 
used to reconstruct the jets. 

• The number of charged particles, A'^^, that are perpendicular in azimuth to the 
leading jet. is a measure of the particle multiplicity associated with the un- 
derlying event. There is no underlying event in the CEP signal events. We use the 
definition adopted in , which assigns charged particles to the underlying event 
if they satisfy 

3 < m - 4>j\ < Y and — < |0fc - 0j| < — (11) 

where (pk is the azimuthal angle of the charged particle and is the azimuthal 
angle of the highest transverse energy jet. 



5.5 Experimental cuts on the final state 

Having defined the experimental variables, we now turn to developing a set of cuts that 
will maximise the significance of the CEP signal. We begin with the variables, Rj and 
Ay. In figure El the distributions for signal events are compared with the hb backgrounds 
from double pomeron exchange and overlap ([p] [X] [p]) events. A cone radius of 0.7 is used 
and the distributions have been normalised to unity. The overlap background extends to 
large values of Rj because there is no correlation between the mass as measured in the 
420m detectors and the central di-jet mass. The DPE background events have a lower 
average Rj because of the necessary presence of pomeron remnants which will often lie 
outside of the jet cones. We define an exclusive-enriched sample by 0.75 < Rj < 1.1, 
Ay < 0.06. The upper limit of the Rj window is chosen to reject overlap background, 
and should therefore be optimised for instantaneous luminosity. 

We now focus on the rejection of inclusive overlap background, [p] [X] [p], utilising the 
charged track multiplicity variables described in section [531 All the distributions shown 
have the cuts on Rj and Ay applied as described above. Charged particles are required 
to satisfy pT > 0.5 GeV and |?7| < 1.75. To take into account the track reconstruction 
efficiency of the ATLAS inner detector for very low pT tracks, 20% of the generated 
charged tracks outside of the di-jets are disregarded at random. Experimentally, charged 
tracks would only be used if they are associated with the di-jet vertex. The approach we 
adopt is to only consider tracks that fall within a vertex window of width ±5 mm around 
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Figure 8: The Rj and Ay distributions are shown in (a) and (b) respectively for the 
signal, [p][&&][p] and DPE [p bb p] backgrounds . The distributions were reconstructed 
using a cone radius of 0.7 after smearing the particles with detector resolution. 



the di-jet vertex. 5 must be larger than the vertex resolution of the inner detector and 
we choose 

6 = ±2a^ (12) 

where az is the vertex resolution of the most poorly reconstructed particle. The vertex 
resolution of the inner detector may be parametrised [53j by 

1 1 5 

a, = 87 + == (13) 

Pt vsm 

Particles with low transverse momentum or large pseudo-rapidity have poorer vertex 
resolution and hence a larger vertex window is applied. 

Figure shows the probability of an inelastic event falling within ^5 mm of the 
primary vertex. From equation [TSj a particle with pT = 0.5 GeV and rj = 2.5 requires 
a vertex window of 6 = ±7.2 mm and therefore, at lO^'^ cm~^ s~^, 80% of events will 
have at least one inelastic pile-up event inside the vertex window. Any cuts on the 
charged particle multiplicity variables will therefore be relatively ineffective. This is the 
reason why the pseudo-rapidity of particles used in the charged particle multiplicity cuts 
is restricted to < 1.75, which allows for a smaller vertex window of ±2.6 mm for a 
particle with px = 0.5 GeV. 

Figure [10] (a) shows the number of charged particles outside of the di-jet system, 
but associated with the di-jet vertex, for the signal, DPE bb and [p][X][p] bb events 
reconstructed using the cone algorithm with a radius of 0.7. The diffractive events (CEP 
and DPE) have few tracks outside of the di-jet system, whereas the inclusive overlap 
background has many more charged particles. Figure[in](b) shows the transverse charged 
particle distributions, N^. The exclusive events are more strongly peaked at zero in this 
distribution. This is as expected because is designed to measure the underlying 
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Figure 9: The probability of at least one inelastic pile- up event falling within a given 
distance, S, of the primary event vertex. 
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Figure 10: The charged track multiplicity outside of the di-jet system, Nc, is shown 
in figure (a). Figure (b) shows the number of charged particles that are transverse 
to the leading jet as defined by equation [TTl In both cases the particles must satisfy 
Pt > 0.5 GeV and \r] < 1.75. Only 80% of the particles are used to replicate ATLAS 
reconstruction efficiency for low pT tracks. 
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Figure 11: The azimuthal angle between the two leading jets for the signal and [p][6&][p] 
background. 

event activity. Similar results are observed for both charged particle distributions using 
the Kt algorithm with Rkj- = 0.55. 

The final central detector variable that we use to reject inclusive overlap events is 
Ac/), the angle between the two leading jets in the event. Central exclusive events have 
no initial state radiation (ISR) and so the di-jet system is produced with no transverse 
boost. This means that the jets will be back-to-back in azimuth, i.e Acf) ~ vr. Inclusive 
overlap events can have a large transverse boost due to ISR. Figure [11] shows the A(p 
distribution for the signal and [p][X][p] bb events. The distribution is again made after 
the Rj and Ay cuts which remove the majority of events that are not back to back due 
to wide angle, final state radiation. We apply the cut vr — \A(f)\ < 0.15 to further enrich 
the exclusive sample. 

5.5.1 Overlap background rejection using fast timing detectors 

Each of the overlap backgrounds can be reduced significantly by matching the di-jet 
vertex to the event vertex reconstructed from the proton time-of-flight information. The 
proposed 420m and 220m detectors at the LHC each have fast timing components for 
this purpose. As we shall see, the timing resolution is critical for rejecting overlap 
background at the highest LHC luminosities. The aim is to achieve a time-of-flight 
measurement with a resolution of lOps [16]. We take this figure as the benchmark in 
what follows. 

The vertex position for overlap events reconstructed by the forward detectors may 
not coincide in space or time with the di-jet vertex reconstructed in the central detector. 
To calculate the timing rejection factor for [p][X][p] events we select three 'event vertices' 
distributed in longitudinal position according to a gaussian distribution of width 4.45cm 
and time according to a gaussian distribution of width 0.18ns [20j . We define the first 
vertex to be the di-jet vertex, and the second and third to be the soft-single diffractive 
vertices that produce the outgoing protons. The fake vertex position Zf as measured by 
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Figure 12: Time-of-flight rejection factor for each type of overlap background, assuming 
a timing resolution of lOps in the forward detectors. 

the forward detectors is given by 



where the t± are the arrival times of the protons measured in the forward detectors in 
the ± z direction. For overlap events made up of two proton-proton interactions (i.e. 
[p] [pX] and [pp] [X] ) a similar method can be used after selecting two event vertices (the 
first of which is the primary di-jet vertex). An event is kept if the fake vertex lies within 
±4.2 mm of the primary vertex (this is ±20" assuming a timing resolution in the forward 
detectors of lOps). 

We must also take account of the possibility that more than one proton might be 
detected in each forward detector. In this case, if any combination of protons produces 
a fake vertex within the ±4.2 mm window then they will survive the vertex cuts. This 
implies that the rejection factor decreases at higher luminosities where the probability 
for additional protons is larger. Figure [12] shows the rejection factor for each type of 
background as a function of luminosity. 

The rejection factors presented in figure [T2] correspond to a lOps timing resolution, 
which is the base design of FP420. It should be noted that if this timing resolution 
is improved, the rejection factors increase linearly with reduction in time-of-flight mea- 
surement, i.e the rejection factors increase by a factor of 5 for 2ps timing resolution. 
Furthermore, it has been proposed that the central detector could provide a time of 
interaction measurement for the primary vertex [20]. In such a case, the time of arrival 
of each proton can be compared to the time of the interaction and this results in a order 
of magnitude increase in rejection for the [p][X][p] overlap background if the central 
interaction is measured to a lOps accuracy. 




(14) 
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6 Results 



6.1 Standard Model Higgs boson 

We focus initially on the symmetric CEP events, i.e. those for which both protons are 
detected in 420m proton detectors at each side of the central detector. CEP signal and 
di-jet background events are generated by ExHuME, and added to the non-CEP di-jet 
backgrounds generated by HERWIG + JIMMY and POMWIG, with forward protons 
added into the event record in the case of overlap events as described in section HI All 
Monte Carlo events will therefore have two forward protons, one on each side of the 
central di-jet system. The 4-momenta of all outgoing particles are smeared to simulate 
the response of the ATLAS detector, as described in section [5. 3. 2[ Jets are found using 
the midpoint cone algorithm with a cone radius of 0.7 (section l5.3p . The top line of table 
[5] shows the cross sections for each sample, after the following cuts on the Monte Carlo 
events: 

• Both forward protons must lie within the geometric acceptance of the 420m forward 
detectors, as described in section [21 with the detectors positioned 5mm from the 
beam. This corresponds approximately to the kinematic range 0.005 < < 0.018, 
0.004 < ^2 < 0.014 and unrestricted in t. 

• The mass of the central system as measured in the forward detectors must satisfy 
80 GeV < M < 160 GeV. 

• Two jets with Eti > 45 GeV, Et2 > 30 GeV. 

The exclusive candidate sample is defined by the following cuts: 

• The di-jet mass fraction 0.75 < Rj < 1.1. 

• The rapidity of the central system Ay < 0.06. 

• The jets must be back-to-back in azimuth, i.e vr — |A(^| < 0.15. 

• The multiplicity of charged tracks with pT > 0.5 GeV and |r/| < 1.75 that are 
associated with the di-jet vertex, Nc < 3 and < 1. In order to simulate the 
ATLAS detector reconstruction efficiency, only 80% of charged particles are used. 

Table [2] shows the effect of the cuts on the signal and background samples. For 
the overlap and DPE backgrounds, only the dominant bb contributions are shown. The 
contribution from the non 6-jet backgrounds are shown after all analysis cuts in table El 
The final cut in table 2 is a mass window, AM, around the Higgs boson peak; 



where Th is the width of the Higgs boson and ctm is the mass resolution of the forward 
detectors. A 120 GeV SM Higgs boson has a very small width, and the size of the mass 
window is therefore dominated by the resolution of the forward detectors. The size of 




(15) 
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Cross section (fb) 


Cut 


CEP 


DPE 


[p][x][p] 


[p][px] 


[pp][x] 




H 


hh 


99 


bb 


bb 


bb 


bb 


Et, £.1, ^2 


0.124 


1.320 


2.038 


0.633 


3.91x10^ 


Ti 

7.33x10^ 


T 

6.29x10"* 


Rj 


0.119 


1.182 


1.905 


0.218 


4.73x10"^ 


85.2 


7.59x10^ 


Ay 


0.010 


1.036 


1.397 


0.063 


2.16x10^ 


1.38 


3.50x10^ 


A$ 


0.093 


0.996 


1.229 


0.058 


6.66x10^ 


0.77 


1.07x10^ 


Nc, 


0.084 


0.923 


0.932 


0.044 


6.49 


0.45 


1.35 


AM 


0.072 


0.070 


0.084 


0.004 


0.59 


0.03 


0.13 



Table 2: Cross section (fb) after applying each cut for the analysis performed using the 
cone algorithm with radius 0.7. The first cut is applied after requiring that both protons 
are tagged at 420m, the mass measured by the forward detectors is between 80 and 160 
GeV and the transverse energy of the leading jet is greater than 40 GeV. The overlap 
backgrounds are defined at high luminosity (10^^ cm~^ s~^) and TOF vertexing has 
been used. 



the mass window is chosen so that more than 80% of the signal is retained. The mass 
window for events in which both protons are tagged at 420m is ±2.4 GeV. 

The largest loss of signal events is caused by the double 6-tagging efficiency of 0.36 
and the acceptance of the forward detectors. It may be possible to improve the b- 
tagging efficiency for central exclusive events due to the clean nature of the CEP vertex. 
In addition, various strategies have been discussed to increase the acceptance of the 
proposed 420m detectors [TB] . 

The final cross sections for both signal and backgrounds are shown in table [3j Note 
that we do not include the trigger efficiency in this table. The final results for the two 
jet algorithms are similar and we use only the cone algorithm for the remainder of this 
paper. The mass window defined by equation [T^ has been applied so that the size of the 
signal and background in the region of interest can be evaluated. At low luminosity, the 
dominant backgrounds are from central exclusive di-jet production. At high luminosity 
the [p][X][p] background becomes dominant. In table [3] we also show the signal and 
background cross sections for the asymmetric events for a distance of approach of (5mm 
+ 2mm) and (3mm + 1.5mm) of the 420m and 220m detectors from the beams. Because 
of the poorer resolution of the 220m taggers, the mass window is increased to ±4.9 

If the cross section for central exclusive SM Higgs production in the b-jet channel is 
~ 2.5 fb, then after experimental cuts one would expect around 0.19 events / fb~^ if the 
events could be triggered at LI and the detectors could be operated at the minimum 
possible distance from the beam. This corresponds to 5.6 signal events after three years 
of data acquisition at low luminosity with a background of 21 events. At high luminosity, 
one would expect a total of 56 signal events with a background of 1032 events after 3 
years data taking. Approximately 80% of the background events at high luminosity are 

^The achievable mass resolution for asymmetric events is slightly different on each side of the IP due 
to a difference in beam 1 and beam 2 optics. Here we quote the average mass resolution. 
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C420-420 (fb) 


^420-220 (fb), CONE 






Kt 


CONE 


5mm/2mm 


3mm/ 1.5mm 




H ^ hh 


0.071 


0.072 


0.038 


0.115 


FvHnA/TF 

JJj A. 1 J. U. iVl JJj 


hh 
uu 


0.076 


0.070 


0.067 


0.203 




fin 

yy 


u.uuu 




n DQi 


n 978 


r DM WIG 


bb 


U.Uii 


U.UU4 


U.UU4 


U.Uio 




ii 


u.uuuo 


u.uuuz 


n nnno 
u.uuuz 


n nnn7 
u.uuu / 


[p][X][p] (L) 


bb 


0.0029 


0.0037 


0.0032 


0.0097 




n_ 


0.0003 


0.0003 


0.0003 


0.0009 


r 1 r-\7"i r i /tt\ 

[p][X][p] (H) 


bb 


0.46 


0.59 


0.46 


1.41 




n 


0.04 


0.05 


0.04 


0.13 


[PPJiXJ (L) 


bo 


U.UUo 


u.uuy 


u.uuy 


U.UZo 


[PPJl^J (J^j 


hh 
00 


n 1 1 


n 1 


n 1 9 


u.oo 




hh 
uu 


0.003 


0.002 


0.002 


0.006 


[P][PX] (H) 


bb 


0.05 


0.03 


0.03 


0.08 


Total bgrd (L) 




0.17 


0.17 


0.18 


0.54 


Total bgrd (H) 




0.81 


0.96 


0.81 


2.48 



Table 3: The final cross sections for the H ^ bb and relevant backgrounds for cone and 
Kt algorithms. The cross sections are defined after a mass window, given by equation 
[T5| has been applied around the Higgs peak. The overlap backgrounds are defined at 
both low (L) and high (H) luminosity, [pp] [X] and [p] [pX] jj backgrounds are negligible 
and not shown. 

from the overlap backgrounds, and it is these backgrounds that would render the signal 
unobservable with a reasonable significance, even if a trigger strategy could be developed 
to retain a high fraction of the events. With the standard experimental configuration, 
therefore, observation of the SM Higgs in the bb channel would not be possible. 

As discussed above, it might be expected that the asymmetric events could be saved 
with a high efficiency by triggering on 220m detectors at LI and reducing the rate at level 
2 using 420m detector hits and timing information. Improving the timing resolution to 
better than lOps, or adding central detector information into the timing system could sig- 
nificantly improve the rejection of the overlap backgrounds as discussed in section [5. 5. 11 
If these improvements are made, a significance of approximately 2 could be achieved. 
Alternatively and additionally, if improvements in trigger efficiency for symmetric tags, 
b-tagging efficiency and acceptance of the forward detectors could be achieved, or the 
cross section turns out to be at the upper-end of theoretical expectations, it may be 
possible to achieve a 3a observation. Since there may well be no other way of measuring 
the SM Higgs b-quark coupling at the LHC, these possibilities should not be dismissed. 
We conclude, however, that observing the CEP of the SM Higgs in the bb decay channel 
will be extremely challenging. We now move on to consider in more detail a scenario in 
which the CEP cross section is predicted to be ~ 8 times higher than the SM prediction. 
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6.2 MSSM Higgs bosons 



In this section, we perform a similar analysis for Higgs bosons produced for the bench- 
mark m^^^ scenario as described in section [3l The cross section for the CEP of the 
lightest Higgs boson in this scenario is enhanced by a factor of 8 relative to the Standard 
Model. 

After the same cuts as for the SM Higgs analysis, except for applying a wider mass 
window of ±5.2 GeV for the symmetric events and ±6.7 GeV for the asymmetric events 
because of the increased width of the Higgs (3.3 GeV), the signal cross section is 0.54 
fb for the 420m ± 420m case, and 0.28 (0.85) fb for the 420m ± 220m case if the 
detectors are 5mm/2mm (3mm/1.5mm) from the beam. This would lead to over 150 
signal events for the symmetric analysis, and over 200 events for the combined analysis 
for an integrated luminosity of 300 fb~^, again assuming a 100 % trigger efficiency. 
Because of the larger mass window, the continuum background within the window also 
increases by a factor of approximately 2-2.5. 

6.2.1 Potential for observation of MSSM Higgs Bosons and mass determi- 
nation 

In order to determine the significance of the MSSM Higgs signal, we construct an event 
sample consisting of the appropriate number of events after 3 years of data taking from 
the signal and background Monte Carlos after the cuts described in section 16.11 We 
assume several different trigger strategies as described in section EH We present results 
for data acquisition at low and high luminosity, i.e. 30 fb~^ of data at a luminosity of 
10^^ cm~^ s~^ and 300 fb~^ of data at 10^"^ cm~^ s~^. In real experimental conditions 
the data will be taken at many different instantaneous luminosities. 

We perform a fit to the pseudo-data with two curves; one assuming that there is a 
signal in the region 100 < M < 140 GeV and one assuming that there is no signal. The 
significance, 5, is then given by 

S={^X^Y (16) 

where A^^ is the difference of the of the fits for signal and null hypotheses. We fit 
the peak region using a gaussian convoluted with a Lorentzian function; the gaussian 
represents the mass resolution of the detectors, which is a known parameter. The width 
of the Higgs boson is left as a free parameter in the fit. We assume that the shape of the 
background will be known, as it can be measured with high statistics with the forward 
detectors. The normalisation of the background is left as a free parameter in the fit 0. 
We repeat this method for 500 pseudo-data sets to obtain the average significance of the 
fits. 

Examples of three pseudo-data sets and fits are shown in figures [13] and [TH In 
figure [13] we consider 60 fb^^ of data taken at 2 x lO'^^ cm^^ s^^ with 420m ± 420m 
tags only. The error bars correspond to \/]V statistical errors. The LI triggers are 

^As a cross-check, the fits were repeated using a simple quadratic fit to the background, and similar 
results were obtained. 
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M (GeV) 

(a) 

Figure 13: A typical mass fit for 3 years of data taking at 2 x 10^'^ cm~^ (60 fb^^). 
The significance of the fit is 3.5a and uses only events with both protons tagged at 420m. 
The LI trigger strategy and analysis cuts are described in the text. 

J25 + MU6 and the rapidity gap trigger (see section [5T]) . Figure (a) is for the same 
experimental configuration, but for 300 fb~^ of data taken at lO^'' cm~^ s~^. The number 
of events increases due to the non-prescaled MU6 trigger, but the signal-to-background 
ratio decreases because the overlap background increases with instantaneous luminosity. 
Figures [TSlandfHl (a) are constructed under the assumption that the 420m detectors have 
a timing resolution of lOps. If the timing resolution is improved or timing information 
from the central detector can be combined with the forward detector information, then 
the overlap background can be effectively eliminated at 10^^ cm~^ s~^ as described in 
section 15.5.11 Figure [E] (b) shows the high-luminosity data set without the overlap 
background. 

Figure [15] (a) shows the significance of the measurement, after 3 years of running 
at instantaneous luminosities between 10^^ cm^^ s^^ and 10^^ cm~^ s~^ assuming that 
the protons are tagged only at 420m. A significance of above 3 is achievable in any 3 
year period, regardless of the instantaneous luminosity if the high rate jet trigger (J25) 
is incorporated at ATLAS in conjunction with the MU6 transverse momentum muon 
trigger. The significance remains approximately constant because although the overlap 
background increases at high luminosity and the pre-scaled jet trigger J25 retains an 
approximately fixed number of events, the non-prescaled muon trigger MU6 retains an 
increasing number of events, and therefore the total number of signal events increases. 
This effect can also be seen in figures [13] and [H] Figure [16] shows the individual contri- 
butions of the muon and jet triggers. The rapidity gap trigger contributes very little to 
the significance at luminosities of 2 x 10^^ cm~^ s~^ and above. 

Figure [15] (b) shows the effect of completely removing the overlap background (i.e. 
using improved time-of-flight rejection). In this case, the more conservative trigger 
combination of a 10 GeV transverse momentum muon (MUlO) and a pre-scaled jet-rate 
of lOkHz (JIO) is enough to achieve a significance of approximately 3 at high luminosity 
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M (GeV) M (GeV) 

(a) (b) 

Figure 14: Figure (a) shows a mass fit for 3 years of data taking at 10^^ cm~^ s~^, 
for events with both protons tagged at 420m. The significance is 3a. The LI trigger 
and analysis cuts are described in the text. Figure (c) shows the same high-luminosity 
data set after removing the overlap background contribution completely, which would 
be possible with improved timing detectors. The significance is 5a 



in a 3 year period. For J25 and MU6, the significance approaches 5. These figures show 
that with improved timing information, running at 10^^ cm~^ s~^ is optimal, whereas 
for the baseline lOps + lOps timing with no information from the central detector, data 
taking at luminosities in the range 2 — 5 x 10^^ cm^^ s^^ is optimal. 

Figure [17] (a) shows the effect of combining the significance of the analysis using 
420m detectors alone with the asymmetric analysis. We use the configuration chosen 
in section [21 where the detectors at 420m are 5mm from the beam and the detectors 
at 220m are 2mm from the beam. Here we choose the JIO + MUlO trigger, since this 
trigger scenario is difficult for 420m detectors alone. The significance rises from ~ 2.5 
to ~ 3.0 for the nominal detector positions (5mm / 2mm) from the beam. Figure [T7I (a) 
also shows the combined significance increases to between 3.5 and 4.0 if the detectors 
are moved to the closest possible distance from the beam (3mm /1.5mm). If the overlap 
background is eliminated then the combined significance rises to between 3.6 and 5.0, 
depending on the distance of approach to the beam. For the J25 + MU6 trigger (not 
shown), the significance after 3 years at 10^"^ cm~-^ s~^ rises to 3.7(4.5), for detectors at 
5mm / 2mm (3mm / 1.5mm) from the beam. If the overlap background is eliminated, 
the significances are 5.5 (7.1) for the J25 + MU6 trigger. 

A further improvement could be made by including the 220m detectors in the level 
1 trigger decision [55] , as discussed in section 15. li Figure [T7| (b) shows the significance 
of the asymmetric analysis if all of the events are retained by the level 1 trigger on 
the 220m detector. For the standard 5mm + 2mm detector positions, the asymmetric 
analysis achieves a significance of 3 at 2 x 10^^ cm~^ s~^. However, at higher luminosities 
the significance approaches 4. If the detectors are moved to 3mm + 1.5mm, then a 
significance of more than 5.0 is achieved at a luminosity of 2 x 10^^ cm~^ s~^. If the 
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Figure 15: Figure (a) shows the significance of the measurement of the 120 GeV MSSM 
Higgs Boson for the model described in the text, using the combined muon, rapidity gap 
and jet rate triggers for the analysis with both protons detected at 420m. Figure (b) 
shows the effect of completely removing the overlap background on two of the trigger 
strategies. 




Figure 16: Figure (a) shows the significance assuming that a low transverse momentum 
muon trigger is used. Figure (b) shows the significance for a pre-scaled jet rate trigger. 
Both significances correspond to an analysis with both protons tagged at 420m and do 
not include the contribution from asymmetric proton tagging. 
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Figure 17: Figure (a) shows the significance of the muon (MUlO), rapidity gap and 
jet rate (JIO) triggers for the combined symmetric and asymmetric analyses. This is 
the more conservative jet rate trigger. Figure (b) shows the contribution from just the 
asymmetric analyses if all of the events are retained by the level 1 trigger. 



overlap background is removed, the significance at high luminosity approaches 11 after 
3 years of data taking. 



7 Uncertainties on the calculated cross sections 

All the signal and background processes are affected by the uncertainty in the calculation 
of the soft survival factor. Whilst this is expected to vary to some extent from process to 
process, to a reasonable approximation it produces a change in the overall normalisation 
of all cross sections, and is not expected to affect the signal to background ratios quoted 
in this paper (although it will of course affect the significance). 

In the analysis presented here it is not possible to quantify the effect of pile- up on the 
jet-finding and the analysis cuts used to define the exclusive sample. Such effects will 
depend on the performance of the central detector as well as the pile-up conditions. We 
now deal with the uncertainties specific to each of the signal and background calculations. 



7.1 Signal 

As discussed in section [3l the cross section for central exclusive production is expected 
to have at least a factor of two uncertainty. This would also apply to the CEP back- 
grounds, since the uncertainty enters primarily in the calculation of the luminosity term, 
which contains the physics of the colour-singlet gluons, and not in the hard sub-process 
cross section calculation of Higgs, bb or di-gluon production. We note that the CDF 
collaboration have observed exclusive di-photon candidate events at the predicted rate 
[54) . It should be noted, however, that the invariant mass of the di-photon system is 
much lower than 120 GeV, and the uncertainties on the calculation may be significantly 
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larger due to non-perturbative effects. CDF have also presented preliminary results on 
the search for the CEP of di-jets [36]. The data are well described if an exclusive di-jet 
component is present in the data at the level predicted by ExHuME (KMR). In sum- 
mary, therefore, the current KMR predictions are consistent with the Tevatron data, 
but there are uncertainties in extrapolating to LHC energies. 

There are many other parameter choices in the MSSM. For lower values of tan/3, the 
cross section is not as large, but neither is the width. As rriA is increased, the heavier 
scalar boson receives an enhanced CEP cross section and the bb decay channel has a 
large branching ratio up to uih ~ 160 GeV. We do not discuss these parameter choices 
and refer to the results in [32] ■ However, we note that the acceptance of a 220m + 
420m forward detector system is approximately constant up to M ~ 200 GeV and so 
the overall conclusions presented in this paper should remain valid for larger rriA. For 
smaller tan/3, the discovery potential will be reduced. 

7.2 Background 

7.2.1 Double pomeron exchange 

The main source of uncertainty in the DPE processes comes from the diffractive parton 
distribution functions at high values of /3, the momentum fraction of the pomeron carried 
by the struck parton. If the value of /3 is large for both partons entering the hard scatter, 
then the reconstructed Rj value will also be large, and this is the region that contributes 
to the background. Thus, the final cross section depends strongly on the high-/3 parton 
distribution of the pomeron. Figure [TH] shows the effect of changing the diffractive pdfs 
on the Rj distribution of the DPE bb events. HI 2006 Fit A and HI 1994 Fit 5 are shown 
against HI 2006 Fit B, which was chosen as the default fit in this analysis. HI 2006 Fit 
A increases the final cross section by a factor of approximately 6. Even with such an 
increase, however, the DPE background stays well below the CEP backgrounds, and the 
overlap background at high luminosity. If HI 1994 Fit 5 is used there is a large increase 
in cross section in the large Rj region (x70). The reason for this is that the DGLAP 
evolution of this diffractive PDF is frozen above = 75 GeV^ , which means that there 
is an increased probability of large /3 values because there is no radiation between the 
pdf freeze-out and the hard scale. The 2006 pdfs include data up to = 1600 GeV^ 
and are evolved outside of this region, and are therefore likely to give a better estimate 
of the background. HI 2006 Fit B has a similar high-/3 behavior to the MRW diffractive 
pdf fits [HIES]. 

7.2.2 Overlap background 

The overlap background calculation of equation [3] requires knowledge of the fraction of 
events that produce one or two forward protons within the acceptance of the forward 
detectors. Table [T] shows that for a single forward proton measured in the 420m region 
(0.002 < < 0.02), there is a consensus between different models and that the fraction of 
events that have such a proton is predicted to be ~1% at the LHC. For larger values of ^, 
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Figure 18: The Rj distribution of the DPE bb events for three different diffractive PDF 
fits. 

the non-diffractive contribution is not well constrained and this gives an uncertainty on 
the overlap cross section with one proton detected at 420m and one at 220m. However, 
the non-diffractive protons typically have larger values of than the diffractive protons 
and so have a smaller chance of producing a central mass in the ~ 120 GeV range. 
In fact, the final cross section for the asymmetric-tagged [p][X][p] background given in 
table El which contains a non-diffractive component, is only 30% larger than without the 
non-diffractive contribution. If the non-diffractive contribution is uncertain to ~ 50%, 
as indicated in table [U then the asymmetric-tag overlap backgrounds are uncertain to 
~ 15%. 

For the [pp][X] events, the uncertainty is larger. In this analysis, we use the total 
DPE cross section predicted by PHOJET to give the [pp] event fraction at the LHC. 
However, PHOJET predicts that a large component of the single forward proton spec- 
trum originates from DPE events. A similar result was observed in [19], where it was 
found that the hard-scattering DPE cross section predicted by PHOJET was much larger 
than that predicted by POMWIG. This implies that the [pp][X] cross section could be 
much smaller at the LHC than we assume here. 

The final uncertainty is in the number of charged tracks produced by an inclusive 
di-jet overlap event at the LHC. This will affect the efficiency of the charged track 
multiplicity cuts on the [pp] [X] and [p] [X] [p] backgrounds. The number of charged tracks 
outside of the di-jet system will be dominated by initial state radiation and multi-parton 
interactions. In this analysis we use HERWIG with JIMMY tuned to CDF data (tune 
A in [56]). A second tune in [56] (tune B) is also able to fit the data. Pythia may be 
used with a variety of different tunes [52] , all of which fit the CDF data. The tunes each 
predict a different number of charged tacks in the transverse region when extrapolated to 
LHC energies. TablejUshows the percentage of the inclusive di-jet backgrounds that pass 
the Nc and A'^^ cuts with the different available tunes. Apart from the Pythia DW tune, 
the rejection factors vary by approximately a factor of 2.5. However, it should be noted 
that the ATLAS detector is capable (although with a reduced efficiency) of measuring 
charged particles that have pT < 0.5 GeV and this could increase the efficiency of the 
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Generator 


Tune 


< 1 


< 1 and 






Nc<3 


HERWIG/ JIMMY 


A 


3.2 


1.0 


HERWIG /JIMMY 


B 


4.1 


1.2 


Pythia 


DWT 


9.5 


2.2 


Pythia 


DW 


14.9 


3.7 


Pythia 


ATLAS 


7.3 


2.3 



Table 4: Percentage of inclusive di-jet backgrounds that pass the charged multiplicity 
cuts. 

charged multiplicity cuts. 

8 Conclusions 

We have shown that the central exclusive production of the lightest scalar Higgs boson, 
for the choice of MSSM parameter space described in section O can be observed with a 
significance of at least 3a in the bb decay channel within 3 years of data taking at the 
LHC, if suitable proton tagging detectors are installed around ATLAS and CMS and the 
current predictions for the CEP cross section are correct. We have evaluated the most 
important backgrounds and shown that they can be rejected with high efficiency using 
a set of exclusivity variables. We have also shown that a fraction of b-jet events can be 
retained by the currently forseen Level 1 trigger hardware, if the trigger strategies we 
outline are adopted. Whilst we have only considered a particular choice of parameters 
in detail, the general conclusions should hold for a wide range of scenarios. As a general 
rule, if the CEP cross section for the production of the Higgs boson is greater than ~ 10 
fb, and the Higgs boson decays predominantly to b-quarks, then the analysis presented 
here should apply if the decay width is not too large. The analysis will also apply to 
any new particle that decays predominantly to b-quarks. 

It is worth speculating what the future experimental strategy might be if a Higgs sec- 
tor such as the m™"^ scenario of the MSSM is discovered at the LHC or Tevatron, and the 
CEP channel proves to be observable with the forward detector configurations currently 
proposed. The largest loss of CEP signal events comes from the limited acceptance of the 
proton detectors (between 25% and 40% depending on detector configuration) and the 
LI trigger efficiency, which is at best around 20% at high luminosity for the strategies we 
consider in this paper. If a hardware upgrade of the LI trigger systems of ATLAS and 
CMS were to increase the trigger latency such that the 420m detector signals could be 
included, then a trigger efficiency of close to 100% could be achieved. Figure[l9](a) shows 
a typical mass fit, with 300 fb~^ of data taken at high luminosity using 420m detectors 
alone (i.e. symmetric events), assuming 100% trigger efficiency. The significance of this 
signal is approximately 6.5a. If it is assumed that, in addition, fast timing detector 
improvements can be made as described in section 16.2.11 then the significance rises to 
nearly lOcr, as shown in figure [19] (b). In this case, a measurement could be made within 
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Figure 19: Figure (a) shows a typical mass fit for the analysis at high luminosity if it 
is assumed that a trigger upgrade allows the 420m detectors to be included in the LI 
trigger at ATLAS and CMS. The significance is approximately 6.5cj. Figure (b) shows 
the same data after removing the overlap contribution. The significance is approximately 
10a. 



100 fb~ . Combining the symmetric and asymmetric analyses increases the significance 
such that a 5a measurement could be achieved within 30 fb~^. The mass of the lightest 
scalar Higgs can be measured with an accuracy of better than 1 GeV. For the 777,™"^ 
scenario considered here, the width of the lightest Higgs is not much larger than the 
mass resolution of the proton detectors, and therefore the extraction of the Higgs width 
from the fits is marginal. In other scenarios with widths in excess of ~ 5 GeV, however, 
a direct measurement of the width would also be possible. We remind the reader that 
observation of any particle in the CEP channel would provide a direct measurement of 
the quantum numbers of the particle. Furthermore, if the pseudo-scalar Higgs is close in 
mass to the lightest scalar Higgs, then CEP would provide an unambiguous separation of 
the two states since the pseudo-scalar cannot be produced. Finally, with such a trigger 
strategy in addition to improved fast timing to reduce the overlap backgrounds, the SM 
Higgs may be obervable in the b-jet channel. 
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